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ABSTRACT 

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud 
(LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like 
halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC- 
mass subhalos, and investigate the lower-mass satellites that were associated with these 
subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z = 0 
satellite population of MW/M31 halos were in a surviving LMC-group before falling 
into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for 
groups with higher mass and/or more recent infall times. Groups of satellites disperse 
rapidly in phase space after infall, and their distances and velocities relative to the 
group center become statistically similar to the overall satellite population after 4 — 8 
Gyr. We quantify the likelihood that satellites were associated with an LMG-mass 
group as a function of both distance and velocity relative to the LMC at z = 0. The 
close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies 
to the LMC suggest that 2 — 4 are likely associated with the LMC. Furthermore, if 
several of these dwarfs are genuine members, then the LMC-group probably fell into 
the MW very recently, < 2 Gyr ago. If the connection with the LMG is established 
with follow-up velocity measurements, these “satellites of satellites” represent prime 
candidates to study the affects of group pre-processing on lower mass dwarfs. 

Key words: Galaxy: formation — Galaxy: halo — galaxies: dwarf — galaxies: Mag¬ 
ellanic Glouds 


1 INTRODUCTION 

The A cold dark matter (ACDM) model predicts an abun¬ 
dance of substructure on all (observable) mass scales. Hun¬ 
dreds of subhalos are predicted to surround Milky Way 
(MW) mass halos dKlvoin et al.lll9^ 1. which can be likened 
to a scaled down versio n of the thousands of substructures 
associated with clusters (iMoore et al.lll999h . This trend, pre¬ 
sumably, continues to smaller mass scales, whereby dwarf 
galaxies can also host several substructures. Recent discover¬ 
ies of dwarf-dwa rf accretion llMartmez-Delgado et al.ll2012l : 

I Belokur^ I 2 OI 3 II present tantalizing observational evidence 
for the existence of such “sub-structure of sub-structure”. 

Despite the hierarchical nature of dark matter halos, we 
generally ignore the possibility that some of the MW satel¬ 
lites may have been part of a group of subhalos before they 
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fell into the Galaxy. The relatively unexplored population of 
sub-subhalos or “satellites of satellites” is strongly linked to 
the most massive structures in the MW halo, as these are the 
potential veh icles that dra g ged in several low mass dwarfs. 
For example. IWetzel et al.l ( 2015l l showed that a significant 
fraction (-^ 30 %) of low-mass subhalos (Mstar < 10 ®Mq) 
likely fell into a MW-type host as a satellite of a more mas¬ 
sive subhalo, and > 50% were in a group before infall. The 
most likely culprit in our own Galaxy is the Large Magellanic 
Gloud (LMG). This massive dwarf already has one obvious 
companion, the Small Magellanic Cloud (SMC), but it likely 
had several other companions in the past. 

Numerous works have attempted to connect th e LMC 
to other known dwarfs in the MW. iLvnden-Belll (Il976l l 
first suggested the idea of a “Greater Magellanic Galaxy”, 
and he later postulated the association of s everal of the 
classi c al dwarfs with the Magellan i c com plex (iLvnden-Belll 
1 19821 : iLvnden-Bell fc Lvnden-Belll Il995h . More recently. 
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iD’Onghia fc La^ (l2008l l suggested that seven of the MW 
satellites co uld have been part o f a late infalling LMC group. 
In contrast, I Sales et al.l (2OIII I use an LMC-analog “case- 
study” in a cosmological simulation to show that most of 
the classical dwarfs show little evidence for an association 
with the LMC. However, the authors do prophetically sug¬ 
gest that faint, previously unnoticed MW satellites could be 
lurking in the vicinity of the Clouds. 

The discovery of very low luminosity galaxies 
{L < 10^L(7)) in t h e MW (e.g, IWillman et al.l l2005l : 
iBelokurov et al.l l2006l . l2007l l has, until recently, been re¬ 
stricted to the Sloan Digital Sky Survey (SDSS) footprint, 
as most of the “ultra-faint” dwarf population have been 
discovered using SDSS imaging. However, uncharted terri¬ 
tory beneath declination 5 = —30° has recently been ex¬ 
plored with the first data release of th e Dark Energy Sur¬ 
vey (DES). Two ind ependent groups jBechtol et al.|[2015l : 
IKoposov et al.l I 2 OI 5 I I unveiled eight and nine candidate 
dwarf galaxies correspondingly in the DES data. Curiously, 
these satellites are mostly of the “ultra-faint” variety and 
are in close proximity to t he LMC. 

In I Wetzel et al.l (l2015h we showed that most of the past 
satellites of an LMC-mass dwarf are likely lower mass sub¬ 
halos, likened to the ultra-faints. Thus, the finding of several 
low luminosity dwarfs in close proximity to the LMC could 
potentially confirm a generic prediction of hierarchical struc¬ 
ture formation. In this work, we use cosmological simulations 
to study the satellite populations of LMC-mass dwarfs ac¬ 
creted onto MW/M31-mass halos, in order to understand 
the potential association between the newly discovered DES 
satellites and the LMC in a cosmological context. 

2 NUMERICAL METHODS 

2.1 ELVIS Simulations 

To study the satellite populations of LMC-mass dwarfs, we 
use ELVIS (Exploring the Local Volume in Simulations, 
ICarrison-Kimmel et al.ll2014l L a suite of 48 high-resolution, 
zoom-in simulations of Milky Way/M31 mass halos (iV/vir = 
1 — 3 X 10^^ Mq). Half of the ELVIS halos reside in a paired 
configuration with separations and relative velocities simi¬ 
lar to those of the MW-M31 pair, while the remainder are 
highly isolated halos mass-matched to those in the pairs. 

Within the high-resolution, zoom-in volumes (span¬ 
ning 2-5 Mpc in size), the particle mass is 1.9 x 10® Mq 
and the Plummer-equivalent force softening is 140 pc. 
ELVIS adopts a cosmological model based on Wilkinso n 
Microwave Anisotropy Probe WMAP7 (iLarson et al.l[20llh . 
with the following ACDM parameters: as, = 0.801, Dm = 
0.266, Da = 0.734, = 0.963 and h = 0.71. See 

iGarrison-Kimmel et al.l (l2014l l for more details on ELVIS. 


2.2 Einding and tracking subhalos 

Dark matter subhalos are identified in ELVIS using the six - 
dimensional halo hnder rockstar (iBehroozi et al.ll2013al L 
and merger trees are constructed usin g the CONSISTENT- 
TREES algorithm llBehroozi et al.l[2013bfl . 

For each subhalo, we assign its primary progenitor 



Figure 1. Distribution of peak subhalo mass, Mpeak (solid gray), 
and stellar mass, Mstar (hashed green), for the 25 satellites with 
Mpeak > 10^^ Mq (masses near that expected for the LMC) in 
the MW/M31 hosts at z = 0 in the ELVIS simulation suite. These 
are the (surviving) “LMCs” that we will discuss in the remainder 
of the paper. The dotted line indicates the observed stellar mass 
of the LMC. 

(main branch) as the progenitor that contains the largest to¬ 
tal mass summed from the subhalo masses over all preceding 
snapshots in that branch. We then compute the maximum 
(peak) mass, Mpeak, ever reached by the main branch of a 
progenitor. 

Throughout this work, we only consider subhalos 
with Mneak > (o r Mstar 2) 5 X 10^ Mq); 

iGarrison-Kimmel et al.l (l2014h found that subhalos down to 
this mass-threshold do not suffer from resolution and nu¬ 
merical disruption issues. 


2.3 Sample of LMC-mass satellites 

We select a sample of LMG-mass satellites of MW7M31 
hosts at z = 0 using all 48 (paired and isolatedjj ha¬ 
los in the ELVIS simulation suite. We select z = 0 satel¬ 
lites with Mpeak > IO^^Mq (Mstar ^ 3 X 10® Mq). This 
lower mass cut is approximately a factor of two lower than 
the LMG mass (Mp eak 2 x 10^^ for Mstar « 2 x 10®, 
Ivan der Marel et al.ll200^ . Stellar masses are estimated for 
the dark matter subhalos using th e Msta r — Mpeak relation 
derived in IGarrison-Kimmel et al.l ll2014h . We exclude the 
satellite in the Sonny & Cher paired simulation that has a 
mass comparable to Its host halo (Mpeak ~ 7 x 10^^ Mq). 

A significant number of the host halos in ELVIS (~ 
50%) do not have any satellites more massive than Mpeak > 
10^^Mq, while some halos have more than one LMC satel¬ 
lite. It is unlikely that relatively low mass MW/M31 ha- 
los with Myir » 10^® Mp) h ost very massive satellites (e.g. 
iBovlan-Kolchin et al.l 120101 ) , so we are biased towards the 
more massive host halos in the ELVIS suite (typically 
(Mvir) ~ 2 X 10 ^®Mq). We also note that some of the 


^ We find no significant differences in our results when just the 
paired halos are used. 
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Figure 2. The fraction (left panel) and number (right panel) 
of satellites of MW/M31 hosts at 2 = 0 that were satellites of 
the surviving LMC dwarfs prior to its infall onto the MW/M31 
host. The peak mass (stellar mass) of the LMCs is shown on the 
bottom (top) x-axis. The color scheme indicates the time since 
infall of the accretion events. Recent and/or massive accretion 
events contribute significant numbers of “satellites of satellites” 
to the 2 = 0 satellite population. 


paired halos were selected to have a satellite compan- 
ion with mass similar to th e LMC (Mstar ^ see 

iGarrison-Kimmel et al.ll2014l) . Thus, our sample is not an 
unbiased (random) selection of MW/M31 mass hosts. 

Our final sample comprises 25 LMC-mass dwarf satel¬ 
lites at « = 0. Fig. [T] shows the distribution of their stellar 
and peak dark matter masses. 

2.4 Finding the (snrviving) satellites of LMC 

dwarfs prior to infall onto the MW/M31 host 

We trace back all z = 0 satellites of MW/M31 hosts and 
identify those that were satellites of a surviving LMC dwarf 
anytime before infall onto the MW/M31 hosts. We assume 
all subhalos with Mp^ak. > W^Mq host luminous galaxies. 
We impose that a subhalo must remain a satellite for at least 
two consecutive time-steps (AT « 400 Myr) in the simn- 
lations to avoid counting particularly transient (and likely 
non-meaningfnl) crossings just within Rvii- 

Note that the LMC dwarfs themselves are now also 
satellites of MW/M31 hosts, but prior to infall are the group 
centrals. 

In total, we identify N = 734 MW/M31 satellites to¬ 
day that were once satellites of LMC dwarfs, where these 
“LMCs” are still intact today. These “satellites of satellites” 
comprise approximately 7% of the surviving satellite popu¬ 
lation of MW/M31 hosts at z = 0, and have typical stellar 
masses of Mstar = 10® — 10® Mq (comparable to the ultra- 
fai nt dwarf galaxy po pulation). This fraction is lower than 
in IWetzel et al.l ll2015ll because we only consider the subset 
of satellites that were satellites of a survivinf^ LMC satellite 
before infall. In this work, we only consider subhalos within 

^ IWetzel et al.l ll2015l l show that approximately half of the overall 
population of group centrals have merged/disrupted by 2 = 0. 


the MW/M31 hosts today, and do not include “field” sub¬ 
halos (i.e. outside of Tvir today) that could have been as¬ 
sociated with an LMC-mass dwarf in the past. It is worth 
noting, however, that these associations do exist, and this 
could be an interesting population to study in future work. 

3 RESULTS 

3.1 Satellites of LMC-mass dwarfs 

Fig. [5] shows the fraction (left panel) and number (right 
panel) of 2 = 0 satellites that were associated with a sur¬ 
viving LMC dwarf before infall onto the MW/M31, against 
the mass of the group central. The color scheme indicates 
the infall timf0 of the accretion events onto the MW/M31 
hosts (blue=recent infall, red=early infall). 

Unsurprisingly, more massive dwarfs have more abun¬ 
dant satellite populations. There is also a dependence on 
infall time onto the MW/M31 host. At a given mass, groups 
accreted more recently have more surviving members at 
2 = 0. 

3.2 Phase-space associations at 2 = 0 

We now consider the current (2 = 0) association in phase- 
space between the LMC dwarfs and their former satellite 
population. In Fig. |3] we show the median velocity (left 
panel) and 3D distance (right panels) between the “LMCs” 
and their past group members as a function of infall time 
onto the MW/M31 hosts. We note that here, and through¬ 
out this work, we use infall time as an orbital constraint in 
our analysis. We find that infall time, rather than either the 
current distance from the MW/M31 host or orbital pericen¬ 
tric distance, shows the strongest relation with the dispersal 
of groups in phase-space (see below and Fig. |3]|. 

After infall, groups become more di spersed in phase- 
space over time (see also lSales et al.l[201lh . For comparison, 
we show the typical average velocity/distance difference be¬ 
tween all satellites of MW/M31 hosts at 2 = 0 and the group 
centrals with the dotted lines. Groups accreted more than 
~ 5 — 6 Gyr ago are well mixed in phase-space today. For 
illustration, the far-right panel of Fig. [3] shows the distribu¬ 
tion of AR for one LMC-group with low median AR. Note 
that this group also has similar dynamical properties to the 
observed LMG-system (see Section [3.311 . 

In the middle right panel, we show the median differ¬ 
ence in configuration space for satellites with AR < 130 
kpc from the group central. This is a rough estimate for the 
ma ximum AR probed by the DES survey around the LMG 
fsee lKoposov et ^l2015l Fig. 20). The proximity of the DES 
satellites to the LMC is striking, especially compared to the 
general population of group members in the simulations. 
This proximity in configuration space not only suggests a 
likely association between the DES dwarfs and the LMC, 
but, if several of these dwarfs are genuine group members, 
then it implies a very recent infall time for the LMC-group. 
Note that the most recent observational constraints on the 
orbits of the LMC/SMC suggest a recent infall time for this 

® Throughout we use “infall time” to define the time since infall 
of a subhalo onto a host halo. 
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Figure 3. The median differences in 3D velocity (left panel) and 3D distance (right panels) at z = 0 between the LMC dwarfs and 
the satellites that were associated with them in a group before falling into the MW/M31 hosts. LMCs that fell in at early times have 
the largest differences in phase space with their satellites today. The colors indicate the number of surviving group members. The black 
dotted lines indicate the average velocity/radial difference between all satellites in MW/M31 hosts and the LMC satellite at 2 : = 0. We 
also show the approximate virial radius for an LMC-mass subhalo with the short-dashed line. The middle right panel shows the median 
difference in configuration space for satellites with AR < 130 kpc. This is a rough estimate for the maximum AR probed by the DES 
survey around the LMC. The colors indicate the fraction of surviving group members that have AR < 130 kpc. The black dashed line 
shows the median distance between the DES dwarfs and the LMC. The furthest right panel shows the distribution of AR for one massive 
group (indicated by the star symbol) with low median AR. 


group (see e. g, Besla et al.||2007|: Bov lan-Kolchin et al.|[201ll : 


iRocha et al. II 2 OI 2 I: iKallivavalil et al.ll2013ll. 

Qualitatively, a pict ure similar t o the above has been 
painted by the study of ISales et alj ll201lll . However, here 
we present the first quantitative evidence of a pronounced 
correlation between z = 0 scatter in the phase-space ex¬ 
hibited by the “satellites of the satellites” for a statistically 
significant sample of accretion configurations. 

In Fig. |4]we show the median velocity and 3D distance 
at z = 0 between the “LMCs” and the MW/M31 satellites 
that were associated with them prior to infall, as a func¬ 
tion of the pericenter distance (top panels) and the z = 0 
distance (bottom panels) from the MW/M31 hosts. Groups 
with smaller pericenters show a larger dispersal in phase- 
space, but the scatter between systems is large, particularly 
Rperi < 100 kpc. For pericenters similar to the observed 
closest approach of the LMC (~ 50 kpc), groups can be 
highly dispersed or can remain in close proximity at z = 0 
(e.g. Median AR ~ 50 — 230 kpc). We also note that LMC- 
mass dwarfs with pericenters close to the observed value of 
the LMC (~ 50 kpc), can have a wide range of infall times 
(^1 — 9 Gyr). The trends with present day distance from 
the MW/M31 hosts are also relatively weak, particularly 
for dwarfs with Dlmc-mw ^ 150 kpc . Finally, we also find 
that the dispersal in phase-space of groups, at least for our 
sample, is only weakly dependent on the orbital eccentricity 
of the LMC-mass dwarfs. 


3.3 Dynamical LMC-analogues 

Our sample of LMC-mass dwarfs is based purely on peak 
subhalo mass (see Fig. [TJ. We now select a subsample of 
these dwarfs based on observational dynamical constraints. 
Figure [5] shows the radial and tangential motion of our sam- 


MlmC 
[Mq X 10“] 

Rperi 

[kpc] 

VR(i?peri) 

[km s“l] 

VT(Rperi) 
[km S“^] 

Tinfall 

[Gyr] 

NsoS 

1.1 

78 

78 

327 

2.6 

41 

1.0 

61 

45 

313 

7.3 

17 

1.8 

54 

106 

341 

1.4 

48 


Table 1. Sample of 3 LMC-mass dwarfs with dynamical prop¬ 
erties similar to the observed LMC system (see main text tor 
details). 


pie of 25 LMC-mass dwarfs at pericenter, where the points 
are colored by pericentr ic distance. The e s timat ed motion 
of the LMC derived by IKallivavalil et ahl ll2013ll is shown 
with the star symbol on this plot (Vr,lmc = 64 ± 7 km s“^, 
Vt.lmc = 314 ± 24 km s“^, Rpeii = 50 kpc). The dashed 
green box highlights 3 dwarfs with similar velocity and peri¬ 
centric distance to that observed for the LMC, and the prop¬ 
erties of these dwarfs are listed in Table [T] Note that one of 
these dwarfs was highlighted in Fig. [3] and has a very low 
median AR. 

In the following subsection, we investigate the probabil¬ 
ity of group membership based on the phase-space dispersal 
of group members at z = 0. Here, we use both the full sam¬ 
ple of 25 LMC-mass dwarfs and the subsample of 3 dwarfs 
that we label as “dynamical analogues”. 


3.4 Likelihood of group-membership 

Fig. 1^ presents the probability of a past association with 
an LMC-mass dwarf as a function of distance (left panel) 
and velocity (middle and right panels) from the massive 
group central. This now includes “interloping” satellites near 
the LMC at z = 0 that were not satellites of the LMC- 
mass host prior to MW/M31 infall. Dwarfs more closely 
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Figure 4. The median differences in 3D velocity (left panels) 
and 3D distance (right panels) at 2 = 0 between the LMC dwarfs 
and the satellites that were associated with them in a group be¬ 
fore falling into the MW/M31 hosts, as a function of the distance 
of closest approach to the MW/M31 hosts (-Rperi? top panels) 
and the z = 0 distance between the LMC-mass dwarf and the 
MW/M31 (DlmC—MW? bottom panels). LMCs with larger peri- 
centers tend to have smaller differences in phase space with their 
satellites today. However, there is a lot of scatter in this rela¬ 
tion, particularly for dwarfs with i?peri ^ 100 kpc. For example, 
dwarfs with 40 < i?peri/kpc < 60 (the approximate pericenter 
of the LMC, illustrated with the gray bands) can have a wide 
range of infall times and phase-space differences. Note that two 
dwarfs with i?peri > 150 kpc are shown as lower limits in these 
plots. The trends with present day distance (DlmC— MW? bot¬ 
tom panels) are weaker. The colors indicate the infall times of 
the LMC-mass dwarfs. 


related in configuration or velocity space are more likely 
to have been group members before infall. For example, 
> 25% of dwarfs within 50 kpc of an LMC-mass dwarf 
were likely associated with this dwarf before infall. We 
show the 3D velocity difference and radial velocity difference 
{AVr = |Vh,SoS — Vr,lmc|) between group members in the 
middle and right-hand panels, respectively. Although, 3D ve¬ 
locity information gives a cleaner distinction between mem¬ 
bers and non-members, radial velocities can also be used to 
assign membership probabilities. 

The significance of infall time onto the MW/M31 host 
is further illustrated in Fig. [G] The red long-dashed and blue 
dotted lines show the fraction of past members as a function 
of radial and velocity difference for late (Tinfaii < 2 Gyr) and 
early (Tinfaii > 5 Gyr) accretion events, respectively. Groups 
accreted a long time ago are now phase-mixed, whereas the 
probability of being associated with a recently accreted LMG 
dwarf is strongly related to the proximity in phase-space. 
We also show the probabilities when we only consider the 
dynamical analogues of the observed LMC with the green 
dashed line. These relations are very similar to the late ac¬ 
cretion event subsample, likely because 2 (out of 3) of the 
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Figure 5. The tangential and radial velocity of our sample of 
LMC-mass dwarfs at pericenter. The points are colored by peri¬ 
centric distance. The square symbol (and err or bar) indicates 
the o bserved velocity components of the LMC llKallivavalil et ^ 
[20 i 3). The dashed-green box indicates 3 dwarfs in our sample 
with similar dynamical properties to the observed LMC. The star 
symbol indicates the dwarf highlighted in Fig. [3] with a very low 
median AR. 


dwarfs in this sample were accreted very recently (see Ta¬ 
ble [TJ ■ Also, it is interesting that the average behavior of 
all recently-accreted LMC-analogues does not strongly dif¬ 
fer from that of the dynamical analogues included here. 

Combining both position and velocity information al¬ 
lows a much easier distinction between previous members 
and the general satellite population. In Fig. [7] we show the 
joint probability distributions based on position and veloc¬ 
ity. The top panels use the full sample of LMC-mass dwarfs, 
the middle panels only include late accretion events, and the 
bottom panels are for the dynamical LMC-dwarf analogues. 
The left-hand panels show that > 90% of dwarfs within 50 
kpc and 50 km s“^ of a LMC-mass dwarf were likely once 
group members. For comp arison, AR — 23 A 2 kpc and 
AVsd = 128 ± 32 km s“^ dKallivavalil et aT]|2013lj for the 
LMC-SMC pair. 

We ensure that the results shown in Fig. [6] and Fig. Q 
are not significantly affected by the “LMCs” in our simula¬ 
tions that have very large pericenters. For example, we find 
that our results are unchanged if we restrict our sample to 
LMC-mass dwarfs with Rperi < 100 kpc. Finally, one may 
expect that the probability of associations between satellites 
is affected by the Galactocentric distance of the LMC-mass 
dwarfs today, as there could be more “interlopers” at smaller 
radii from the MW/M31. However, we find no significant 
differences if we only include dwarfs inside (or outside) of 
D — 100 kpc from the host center. 

We can use these relations shown in Fig. and Fig. [7] 
to estimate the probability that the DES candidate dwarfs 
were once satellites of the LMC. The estimated probabilities 
are listed in Table [51 We also give the sum of these proba¬ 
bilities, which provides a rough estimate of the number of 
these dwarfs that are “satellites of satellites”. Using only 3D 
coordinate information, we find that two of the DES dwarfs 
were once satellites of the LMC. If we assume that the LMC- 
group fell in very recently (Tinfaii < 2 Gyr), or restrict our 
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Figure 6. The fraction of all MW/M31 satellites at ^ = 0 that were a satellite of a surviving LMC-mass dwarf before infall onto the 
MW/M31 host as a function of 3D distance (left panel) and velocity (middle and right panels) difference from the LMC today. The 
long-dashed red and dotted blue lines are for groups accreted recently (Tinf^n < 2 Gyr) and early (Tjnfaii > 5 Gyr), respectively. Only 
groups accreted recently show a close-proximity in phase-space at 2 = 0. The dashed-green lines are for the dynamical LMC-analogues 
(see Section 13.311 . 


sample to dynamical analogues of the LMC, then this num¬ 
ber rises to four. Fig. [7] show that the inclusion of velocity 
information will enable a clearer distinction between mem¬ 
bers andnion-unernbers in the future. _ _ 

I Simon et al.l (l2015l ) and IWalker et al.l ll2015l ) recently 
spectroscopically confirmed that the Ret 2 dwarf candidate 
is indeed a dwarf galaxy. Although this dwarf is in close 
proximity to the L MC (AR = 23.9 k pc), t he radial ve¬ 
locity measured by I Simon et al.l ll2015l) and IWalker et al.l 
(120151) for this dwarf is more disparate (|AVr| = 160 km 
s“^Jj. This dwarf is shown by a purple star in Fig. [7] With 
distance information alone we estimated Plmc sat = 0.38, 
Plmc sat (Tinfaii < 2Gyr) = 0.65 and Plmc sat (dyn. ana¬ 
logues) = 0.77, but this additional radial velocity informa¬ 
tion changes the probability of once being a satellite of the 
LMC to Plmc sat = 0.28, Plmc sat(rinfaii < 2Gyr) = 0.57, 
Plmc sat (dyn. analogues) = 0.84 respectively. Finally, we 
also includ e the recent measurem ent of the radial velocity 
of Hor 1 bv iKoposov et aP (l2015l) . to estimate probabilities 
of Plmc sat = 0.36, Plmc sat (Tinfaii < 2Gyr) = 0.74 and 
Plmc sat (dyn. analogues) = 0.77. 

Note that we have not taken into account the presence 
of the SMG in the above analysis, and it is worth point¬ 
ing out this potential caveat. It is beyond the scope of this 
work to quantify the affect of the LMG-SMG interaction on 
the orbits of lower-mass LMC satellites, but this could be a 
worthwhile avenue for further study. 


4 CONCLUSIONS 

We used the ELVIS simulation suite to study the sur¬ 
viving satellite population of LMC-mass dwarfs accreted 
onto MW/M31 mass halos. A sample of 25 LMC-mass 
(ALpeak > IO^^Mq) 2 = 0 Satellites of MW/M31 hosts are 
selected, and we find the lower mass dwarfs that were as¬ 
sociated with these massive dwarfs before they fell into the 

^ Note that this is the difference in line-of-sight velocity between 
Ret 2 and the LMC in the Galactic rest frame. 


Name 

AP 

[kpc] 

-Plmc sat 

-Plmc sat 
(Pnfall < 2 Gyr) 

-Plmc sat 
dyn. analogues 

Reticulum 2 

23.9 

0.38 

0.65 

0.77 

Eridanus 2 

337.4 

0.02 

0.01 

0.04 

Horologium 1 

38.5 

0.31 

0.57 

0.60 

Pictoris 1 

70.0 

0.19 

0.41 

0.29 

Phoenix 2 

54.3 

0.23 

0.49 

0.41 

Indus 1 

80.0 

0.18 

0.37 

0.22 

Crus 1 

92.8 

0.16 

0.31 

0.13 

Eridanus 3 

48.2 

0.26 

0.52 

0.48 

Tucana 2 

36.7 

0.32 

0.58 

0.62 


Total: 

2.0 

3.9 

3.6 


Table 2. The nine candidate dwarf galaxies from lKoposov et al.l 
(120151) . We give the dwarf name, 3D distance from the LMC, and 
estimated probability of once being a satellite of the LMC based 
on this distance. The probabilities are computed using all LMC- 
mass satellites (Plmc sat)i LMC-mass satellites with recent infall 
times (Plmc sat[Tinfall < 2 Gyr]) and LMC-mass satellites with 
similar dynamical constraints to the observed LMC (Plmc sat 
[dyn. analogues]). 

MW/M31 hosts. Our selection is motivated by the recent 
discovery of nine candidate dwarf galaxies in the vicinity of 
the LMC/SMC group. We also identify a subsample of 3 
LMC-mass dwarfs with similar dynamical properties to the 
observed LMC-system, and we compare these with the over¬ 
all sample. Our main conclusions are summarized as follows: 

• Recent, massive accretion events likely “dragged in” a 
significant number of MW/M31 dwarfs. Typically, 7% of the 
surviving 2 = 0 satellite population were once associated 
with surviving LMC-mass dwarfs, but this fraction can vary 
between 1% and 25% depending on the mass and infall time 
of the group central. 

• Groups of dwarfs quickly disperse in phase-space af¬ 
ter infall onto MW/M31 mass hosts. We find that 2 = 0 
MW/M31 satellites that were once satellites of a surviving 
LMC dwarf can typically have large differences in velocity 
or configuration space relative to their group central if they 
fell into the MW/M31 host more than 5 Gyr ago. 

• The proximity of the candidate DES dwarfs to the LMC 
suggests that: (1) several were likely satellites of the LMC at 
some point in the past, and; (2) if they are genuine “satellites 


© 2015 RAS, MNRAS 000, HHz] 














































Satellites of LMC-Mass Dwarfs 7 




0 50 100 150 200 250 


AR [kpo] 



0 50 100 150 200 250 

AR [kpo] 

250 

200 

150 

100 

50 

0 



\ . \. 

I I 


1! 

\ 

I 

\ 

■O' \ 

I 


x\ 

\ ^ 


1^1 \ 

\ ° 


I \ 

\ \ 



\ \ 


w 

E 

> 

<I 


0 50 100 150 200 250 

AR [kpo] 


0 50 100 150 200 250 

AR [kpo] 



AR [kpo] AR [kpo] 


Figure 7. The fraction of all MW/M31 satellites at 2 = 0 that 
were a satellite of a surviving LMC-mass dwarf before infall onto 
the MW/M31 host as a function of 3D distance and velocity dif¬ 
ference from the LMC. The contours levels give fractions of 0.1, 
0.3, 0.5, 0.7 and 0.9 with increasing line thickness. The gray bands 
indicate the range of AR for eight candidate DES dwarfs (exclud¬ 
ing Erl 2). The purple stars ind icate the (spectroscopically con¬ 
firme d) Ret 2 and Ho r 1 dwarfs llSimon et al.ll2015l : IWalker et al.l 
I2OI5I : IKodosov et al.l[2015^ .All groups are shown in the top pan¬ 
els, recently accreted groups (Tinfaii < 2 Gyr) are shown in the 
middle panels and the dynamical LMC-analogues are shown in 
the bottom panels. 
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of satellites”, then the LMC-group was likely accreted very 
recently (< 2 Gyr) for these dwarfs to retain such a close 
proximity in configuration space with the LMC. Distance 
information alone suggests that two to four of the newly dis¬ 
covered DES dwarfs were satellites of the LMC-group before 
infall. 

• The DES dwarfs that were/are satellites of the LMC 
could be prime candidates to study the affects of group pre¬ 
processing. If the LMC-group fell in very recently onto the 
MW, then the members may have spent a significant amount 
of time in this group before joining the MW. In future work, 
we plan to study the affects of group pre-processing in more 
detail. 
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